Atomic Quantum Simulation of Dynamical Gauge Fields coupled to Fermionic Matter: 
From String Breaking to Evolution after a Quench 
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Using a Fermi-Bose mixture of ultra-cold atoms in an optical lattice, we construct a quantum 
simulator for a U(l) gauge theory coupled to fermionic matter. The construction is based on 
quantum links which realize continuous gauge symmetry with discrete quantum variables. At low 
energies, quantum link models with staggered fermions emerge from a Hubbard-type model which 
can be quantum simulated. This allows us to investigate string breaking as well as the real-time 
evolution after a quench in gauge theories, which are inaccessible to classical simulation methods. 



Recently, the condensed matter and atomic physics 
communities have mutually benefited from synergies 
emerging from the quantum simulation of strongly cor- 
related systems using atomic setups HUH- In particu- 
lar, physically interesting quantum many-body systems, 
which can not be solved with classical simulation meth- 
ods, are becoming accessible to analog or digital quantum 
simulation with cold atoms, molecules, and ions. In the 
future, quantum simulators may also enable us to ad- 
dress currently unsolvable problems in particle physics, 
including the real-time evolution of the hot quark-gluon 
plasma emerging from a heavy-ion collision or the deep 
interior of neutron stars [5]. 

The challenge on the atomic physics side is to find 
a physical implementation of gauge theories with cold 
atoms, and to identify possible atomic setups represent- 
ing dynamical gauge fields coupled to fermionic matter. 
Below we provide a toolbox for a U(l) lattice gauge the- 
ory (LGT) using atoms in optical lattices [TJ [3]. Here 
fermionic atoms represent matter fields. They hop be- 
tween lattice sites and interact with dynamical gauge 
fields on the links embodied by bosonic atoms. The LGT 
to be implemented is a so-called quantum link model 
(QLM) [BHS] , where the fundamental gauge variables are 
represented by quantum spins. QLMs extend the con- 
cept of Wilson's LGT [9 J. In particle physics they provide 
an alternative non-perturbative formulation of dynamical 
Abelian and non- Abelian gauge field theories [8j [TUJ [11] . 
QLMs are also relevant in condensed matter contexts, like 
spin liquids and frustrated systems [T2HI4] . Their Hamil- 
tonian formulation provides a natural starting point for 
quantum simulation protocols based on atomic gases in 
optical lattices jT5j[T9l . We will illustrate atomic quan- 
tum simulation of an Abelian QLM in a ID setup, demon- 
strating both dynamical string breaking and the real- 
time evolution after a quench, which are also relevant in 
QCD. The quantum simulator discussed below makes the 
corresponding real-time dynamics, which is exponentially 
hard for classical simulations based on Wilson's paradigm 



[20] . accessible to atomic experiments. 

Cold quantum gases provide a unique experimental 
platform to study many-body dynamics of isolated quan- 
tum systems. In particular, cold atoms in optical lat- 
tices realize Hubbard dynamics for both bosonic and 
fermionic particles, where the single particle and inter- 
action terms can be engineered by external fields. The 
remarkable experimental progress is documented by the 
quantitative determination of phase diagrams in strongly 
interacting regimes, the study of quantum phase tran- 
sitions, and non-equilibrium quench dynamics [2TH25] . 
One of the most exciting recent developments are syn- 
thetic gauge fields with atoms, which promises the re- 
alization of strongly correlated many-body phases, such 
as, e.g., the fractional quantum Hall effect with atoms 
[26H31] . A fermion that is annihilated by ip y and recre- 
ated by ip x at a neighboring site x, which propagates 
in the background of a classical Abelian vector poten- 
tial A gives rise to the hopping term ip x u xy ip y with 
u xy — ex.-p(icp xy ). Hopping between the adjacent lat- 
tice sites x and y accumulates the phase ip xy = j x dl - A. 
The hopping term is invariant against U(l) gauge trans- 
formations A! = A — Va[32 , 33 . When a fermion hops 
around a lattice plaquette (wxyz), it picks up a gauge 
invariant magnetic flux phase exp(z<£) = u wx u xy u yz u zwi 
with = f d 2 f • V x A. We emphasize that these syn- 
thetic gauge fields are c-numbers mimicking an external 
magnetic field for the (neutral) atoms. 

Instead, here we are interested in dynamical gauge 
fields as they arise in particle physics [34]. The corre- 
sponding fundamental bosonic degrees of freedom U xy 
are no longer related to an underlying classical back- 
ground field A, but represent quantum operators as- 
sociated with the lattice links. The hopping of the 
fermions is now mediated by the bosonic gauge field 
via the term iplU xy ip yi which is invariant under local 
changes of matter and gauge degrees of freedom U' xy = 
VW xy V = ex.p(ia x )U xy ex.p(-ia y ), ijj r x = V^ X V = 
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exp(ia x )i/; x , V H x exp (ia x G x ), and G x = ^\^ x - 
E< - • Here E x,x+i is an electric field op- 

erator associated with the link connecting x and y = x+i, 
where i is a unit-vector in the i-direction. G x is the 
generator of gauge transformations (see [33] for a de- 
tailed discussion). Gauge invariant physical states must 
obey Gauss' law, G x \^/} = 0, which is the lattice vari- 
ant of V • E = p = i/j^i/j. To ensure gauge covariance 
of U xy , it must obey [E xy ,U xy ] = U xy . The Hamil- 
tonian representing the electric and magnetic field en- 
ergy of a compact U(l) LGT, H = ^E( xy ) E ly ~ 
4^ T,( WX yz) ( U wxU xy U yz U zw + h.c), is gauge invariant, 
i.e. [H, G x ] = 0. In Wilson's LGT, the link variables 
U xy = ex.-p(icp xy ) G U(l) are still complex phases, and 
E xy = —id/dip X y. Since U xy is a continuous variable, 
which implies an infinite-dimensional Hilbert space per 
link, it is not clear how to implement it in ultra-cold 
matter, where one usually deals with discrete degrees of 
freedom in a finite-dimensional Hilbert space. 

Quantum link models offer an attractive framework 
for the quantum simulation of dynamical gauge fields 
[8|HD|HT]. They extend the concept of a LGT to systems 
of discrete quantum degrees of freedom with only a finite- 
dimensional Hilbert space per link. In contrast to the 
Wilson formulation, QLMs resemble a quantum rather 
than a classical statistical mechanics problem. The re- 
lation [E xy ,U xy ] = U xy is then realized by a quantum 
link operator U xy = S£ y which is a raising operator for 
the electric flux E xy = S xy associated with the link con- 
necting neighboring lattice sites x and y. A local SU (2) 
algebra is generated by a quantum spin S xy with just 
25 + 1 states per link [33]. We will consider quantum 
links with S = \ or 1. In the classical limit S —> oo 
QLMs reduce to the Hamiltonian formulation [35 , 36] of 
Wilson's LGT. 

The implementation of quantum link models in ultra- 
cold matter requires the realization of a gauge invariant 
Hamiltonian accompanied by the corresponding Gauss 
law. Here, we present a general procedure to obtain U(l) 
QLMs including both gauge and matter fields. To illus- 
trate our method, we focus on a simple example, a ID 
U(l) QLM coupled to so-called staggered fermions with 
the Hamiltonian 

H = -t^2 [^l u x,x+i^x+i + h.c] 

X 

+ m X)(-l)^x + yE <x+i- (!) 



Here t is the hopping parameter (see Fig. [IJi), m is the 
fermion mass, and g is the gauge coupling. In this case, 
the gauge generator is given by G x = G x + \ [(— l) x — 1]. 
Staggered fermions are analogous to spinless fermions at 
half- filling in condensed matter physics. The correspond- 
ing vacuum represents a filled Dirac sea of negative en- 
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FIG. 1. [Color online] a) Correlated hop of a fermion as- 
sisted by U x ,x+i = Sxx+i consistent with Gauss' law in a 
QLM with spin S — 1. b) Realization of the process in a) 
with bosonic and fermionic atoms in an optical super-lattice 
(see text), c) Breaking of a string connecting a static QQ pair: 
from an unbroken string (top), via fermion hopping (middle), 
to two mesons separated by vacuum (bottom), d) From a 
parity- invariant staggered flux state (top), via fermion hop- 
ping (middle), to the vacuum with spontaneous parity break- 
ing. 



ergy states. For S = 1, t = 0, and m > the vac- 
uum state has E x ^ x +i = and ^ x ^ x = \ [1 — (— 1)^]. 
The corresponding vacuum energy of a system with L 
sites is Eq = —mL/2. The above Hamiltonian resembles 
the Schwinger model [37] . For S = 1 it shares the non- 
perturbative phenomenon of string breaking by dynam- 
ical qq pair creation with QCD [38]. An external static 
quark-anti-quark pair QQ (with the Gauss law appro- 
priately taken into account) is connected by a confining 
electric flux string (Fig. [TJ^, top), which manifests itself 
by a large value of the electric flux. For t = 0, the energy 
of this state is String — Eq = g 2 (L — l)/2, and the flux is 
given by (%2 X E x ^ x+ i) = — L + 1. At sufficiently large L, 
the string's potential energy is converted into kinetic en- 
ergy by fermion hopping, which amounts to the creation 
of a dynamical quark-anti-quark pair qq. In this process, 
which is known as string breaking, an external static anti- 
quark Q pairs up with a dynamical quark to form a Qq 
meson. For t = 0, the resulting two-meson state of Fig. 
[1J (bottom) has an energy £ mesons — Eq = g 2 + 2m and 
a small flux (%2 x E XjX +i) = —2. The energy difference 
^string — ^mesons = 9 2 (L — 3)/2 — 2m = determines the 
length L = 4m / g 2 + 3 at which the string breaks. 

Another non-perturbative process of interest in parti- 
cle physics is the real-time evolution after a quench. In 
particular, the quark-gluon plasma created in a heavy- ion 
collision quickly returns to the ordinary hadronic vac- 
uum. This is accompanied by the spontaneous break- 
down of the quark's chiral symmetry. The dynamics af- 
ter a quench can be quantum simulated by using the 
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representation for the electric flux (which mim- 



ics the Schwinger model at vacuum angle = tt [37]). 
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FIG. 2. [Color online] Schematic view of the optical super- 
lattices for one fermionic and two bosonic species 1 and 2 
(model I), a) Species 1 can hop between an even site x and the 
odd site x — 1, while species 2 can hop between x and x + 1. 
b) Illustration of various contributions to the Hamiltonian. 
Fermions and two-component bosons have on-site repulsions 
Uif = U2F = U12 = 2U, while bosons of the same species 
have Un = U22 = 2U + g 2 /2 - t 2 B /U. The offsets of the 
bosonic and fermionic super-lattices are 2U± = 2U2 = 2U and 
2Uf = 2(U + m), respectively. If the fermion hops to the 
left, it picks up the energy offset 2U from a boson of species 
2 which simultaneously tunnels to the right. 

In that case, like chiral symmetry in QCD, for m > 
parity is spontaneously broken, at least for small t, for 
more details see [33]. A quenched parity-invariant stag- 
gered flux state, which evolves into the true vacuum with 
spontaneous parity breaking, is schematically illustrated 
m Fig. [TJi. In this case, the electric flux represents an 
order parameter for spontaneous parity breaking, which 
is expected to perform coherent oscillations. This is sim- 
ilar to the time evolution after a quench starting from a 
disoriented chiral condensate in QCD [39] . 

The realization of an atomic LGT simulator requires: 
(i) the identification of physical degrees of freedom to 
represent fermionic particles and bosonic quantum link 
variables; (ii) to impose the Gauss law in order to re- 
move the gauge variant states; and (Hi) to design the 
desired dynamics in the gauge invariant subspace. Below 
we develop a rather general atomic toolbox to implement 
U(l) lattice gauge models coupled to matter fields based 
on mixtures of cold fermionic and bosonic atoms in opti- 
cal lattices. Within this toolbox, we consider two differ- 
ent microscopic realizations in terms of Hubbard models, 
model I and II. Below we present in some detail the con- 
ceptually simpler model I (see Fig. [2|, which assumes 
two-component bosons representing gauge fields. Model 
II, discussed in the [33], assumes one component bosons 
with magnetic or electric dipolar interactions; it offers 
better scalability and experimental feasibility. Our con- 
cepts generalize immediately to experiments in 2D and 
3D, and to fermions with spin [33] . 

(i) The spin S = representing the quantum 



link can be realized with a fixed number N — 2S of 
bosonic atoms in a double well potential with tunnel cou- 
pling (Fig. ^p). An optical super-lattice gJ3 ED (Fig. [2} 
provides an array of double wells with different depths, 
and a Mott insulator phase of bosons allows loading with 
the desired number of atoms N. For two neighboring 
sites x and x + 1, with b x and b x+1 denoting the bo- 
son destruction operators in the corresponding wells, we 
define a Schwinger representation for the quantum link 

U XtX+1 = b? +1 F x , E XtX+1 = \ (Ox+i " K^l) ■ (2) 

The electric flux is related to the population difference of 
the two sites. Here the bosonic species index a = 1, 2 dis- 
tinguishes between links originating from even and odd 
sites x. Eq. Q requires that each boson can tunnel 
only to one specific neighboring site, based on a term 
hxx+i ~ ~ ^b^x+i^ + n - c - The number of bosonic atoms 
is conserved locally on each link. In [33] we discuss 
model II with just a single bosonic species, by encod- 
ing a in the geometric location of the bosons to the left 
or to the right of the site x. We now also add spinless 
fermionic atoms at half-filling to our super-lattice setup, 
which can hop between neighboring sites based on the 
term h x x+1 = — tpV^+iV^c + h.c. (ii) Gauss law: Using 
+ K+1K+1 ~ ^5 ^he gauge generator reduces to 

G x =n F x +nl+nl-2S+ l -[(-lT-l}. (3) 

Here n x counts the atoms of type a — F, 1, 2. Up to an 
x-dependent constant, G x thus counts the total number 
of atoms at the site x. To impose the Gauss law, we 
will consider interaction terms which can be rewritten in 
the form U G x as the dominant term in the Hamiltonian, 
so that all gauge variant states are removed from the 
low-energy sector. This is reminiscent of the repulsive 
Hubbard model for a Mott insulator [1 . In this sense, 
the gauge invariant states (which obey n x +n x + n x = 
25+^ [1 — (— l) x }) can be viewed as "super-Mott" states. 
(Hi) It is well known that, for large on-site repulsion, 
the Hubbard model reduces to the t-J model [42 . We 
now induce the dynamics of a U(l) QLM in a similar 
manner, by considering the ID microscopic Hamiltonian 

H = Ex hlx+i +E X hlx+i+m £x("ir< +U Ex Gl 
Up to an additive constant, it can be expressed as 

H = ~t B h V h l+l ~ f Bj2 b ^ b l+l - l F ^2^x+l 

x odd x even x 

+ h.c. + J2 n a x U af3 4 + J2(-l) X U a n%- (4) 

x,a>,f3 x,a 

The last two terms describe repulsive on-site interactions 
as well as super-lattice offsets, and form the basic build- 
ing block for the Gauss constraint U ^2 X G x . The various 
contributions to the Hamiltonian are illustrated in Fig. 
[2]o. The QLM of Eq. ^ with t = t B t F /U emerges in 
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FIG. 3. [Color online] a) Flux configuration in the ground 
state of EqJZ] compared to the QLM for S — \ obtained by 
exact diagonalization of an L — 8 site system. The parame- 
ters of the QLM (in units of £f = £b = 1) are t — 0.05, 5f — 
—0.05 (see [33 ), and m — —0.2,0,0.2 (squares, crosses, and 
circles). The corresponding microscopic parameters are U = 
20 and m — —0.2,0,0.2 (dashed-dotted, dashed, and solid 
lines), b) accuracy of the effective gauge invariance para- 
meter G — \{G X )\/L in the microscopic realization as a 
function of tp/U.c, d) Real-time evolution of the total electric 
flux E — E x , x +i obtained by exact diagonalization of the 
QLM with L — 16. c) For S = 1 (solid line) string breaking is 
illustrated, starting from the initial state at the top of Fig. [I]:, 
and approaching the corresponding vacuum expectation value 
(dashed-dotted line) of E = J2 X E x , x +i (< 
0, Sjp — — y/2t\ critical breaking length L, 
d) For S = \ we show the evolution after a quench, starting 
from the initial state at the top of Fig. [TJl. The flux order 
parameter performs coherent oscillations whose period and 
strength strongly depends on m (m/t = 0.6(0.9) for dashed 
(thick) line, S F = 10 £). 



2 = y/2t > 0, m = 
= when t 0). 



second order perturbation theory, if one tunes the para- 
meters to the values listed in Fig. [2]o. The offsets U a give 
rise to an alternating super-lattice for both the fermions 
and the bosons. In analogy to super-exchange interac- 
tions [41], energy conservation enforces a correlated hop 
of the fermion with the spin-flip on the link, thus realizing 
the term — tip x U XyX +iil) x +i. This is the key ingredient for 
the coupling of fermions and quantum links. Addition- 
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ally, a gauge invariant term 8p ^2 X 

is also generated [33 . The reduction of the microscopic 
model of Eq. Q to the QLM of Eq. has been verified 
both at the few- and many-body level, is schematically 
illustrated in Fig. [3^i-b and extensively discussed in [33] . 

We have performed exact diagonalizations on small 
system sizes to quantitatively show the physical phe- 
nomena of string breaking and the dynamics after a 
quench which can be observed in an experiment. The 
main results are presented in Fig. [3]>d. For 5 = 1, 
we evolve a string state initially prepared as in Fig. [1J 
under Hamiltonian parameters such that the separation 
between charge and anti-charge is larger than the char- 



acteristic scale for string breaking L = 4m / g 2 + 3. In- 
deed, the large negative electric flux initially stored in 
the string quickly approaches its vacuum value, illustrat- 
ing the string breaking mechanism. For S = |, Fig. [SJi 
also shows the time evolution after a quench, starting 
from the parity- invariant state at the top of Fig. [T]i. 
In fact, the electric flux, which is an order parameter 
for spontaneous parity breaking, displays coherent os- 
cillations, reminiscent of a disoriented chiral condensate 
in QCD [39]. A general experimental implementation, 
which will require three basic steps (preparation of an 
initial gauge invariant state, evolution via quantum link 
dynamics, and measurement of relevant physical observ- 
ables), is discussed in the supplementary materia [33], 

In the present work, we have proposed a quan- 
tum simulator of lattice gauge theories, where bosonic 
gauge fields are coupled to fermionic matter, allowing 
demonstration experiments for phenomena such as time- 
dependent string breaking and the dynamics after a 
quench. While the basic elements behind our model have 
been demonstrated individually in the laboratory, the 
combination of these tools and the extension to higher di- 
mensions remain a challenge to be tackled in future gen- 
erations of optical lattice experiments. While building a 
QCD quantum simulator to address questions related to 
non-zero baryon density and real-time evolution remains 
a long term goal, we see no fundamental obstacles on 
the atomic physics side, but rather a long list of chal- 
lenges such as incorporation of multi-component quark 
fields and non-Abelian plaquette terms in higher dimen- 
sions. A realistic pathway will be the investigation of 
increasingly complex (quantum link) models in an inter- 
play between theory and experiment, with the short term 
goals of extending the present study to higher dimensions 
and in particular non-Abelian gauge field models. 
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SUPPLEMENTARY INFORMATION 



the diagonal Hamiltonian 

2\ 



Symmetries of the U(l) quantum link model H v = fu-\-^\ E «) 2 + 2/7^ n\n\ + U^(-l) x n a x 



In this section, we briefly review the basic symmetry 
properties of the U(l) quantum link model of Eq. (3) of 
the main text. 

1. As in any gauge theory, the Hamiltonian is invari- 
ant against local symmetry transformations. In 
this case, it commutes with the infinitesimal U(l) 
gauge generators (additional details are given in the 
second supplementary material file): 



G x = + ^ [(-!) X - 1] - £x,x+i + (5) 



2. The parity transformation P is implemented as 



U x ,x+1 — C^ljc-i - x , P Ex,x+l — ~E—x—l,—xt (6) 



3. while charge conjugation C acts as 



C U XyX +i = Ul +1 ^ x+2 i C E x ^ x +i = — E x +i^ x +2- (7) 



4. For m — the Hamiltonian also has a Z(2) chi- 
ral symmetry which shifts all fields by one lattice 
spacing, 

x ^x =i>x+i, =^i+i, 

X U x , x +\ = x E x , x +\ = E x+ i jX+ 2- (8) 

However, this symmetry is explicitly broken when 
one imposes the Gauss law G x \^f) — 0. 



Model I: Quantum link model emerging from a 
Hubbard-type model 



In this section, we sketch the main steps to reduce the 
microscopic Hubbard model, Eq. (6) of the main text, 
to an effective quantum link model at low energies using 
second order perturbation theory. We are interested in 
the scenario where the largest energy scale U is given by 



2U n^n x + (U + m)J2(-r 

X, (7=1,2 X 

J2 



x F 



U 



E ( E x-l,x + E l,x+l) 



+ U E (— l) x (E x -i jX — E XiX +i) 

X 

+ 2/7 E [i> x i> x (Ex-i,x — Ex,x+i) — Ex-i,xEx,x+i] 

X 

+ {u+m) j2(-i) x ^ x 

X 

= f/ E^ + Y^^ +1+ra ^ ( " ir ^' (9) 

XX X 

The values g 2 and m are small compared to U > 0, i.e. 
g 2 , \m\ <C U, but they are still relevant in the induced 
quantum link model. 

The term to be generated in second order perturbation 
theory is the correlated hopping of fermions mediated by 
the quantum link (represented by a quantum spin). It 
appears as an effective interaction induced by the previ- 
ous Hamiltonian and the perturbation terms 

AH = -t F (fx+i^x + 1>l*l>x+i) 

X 

-t B [W+i+Oi 



odd 

-t B E \b 2 M +1 + b 2 J +1 bl 



(10) 



To second order in tp and the effective Hamiltonian 
reads 



2 

tpts 



t 2 

<_B_ 

u 



E hfc 

,x+l Vx 



u 



J2^ x (l -V'l+iV'x+i) • 



(ii) 



The last term proportional to 5p = t 2 F /U was not present 
in the original quantum link model Hamiltonian. This is 
no problem, because this term is also gauge invariant, and 
could have been added to the quantum link Hamiltonian 
from the beginning. 

To test the reduction of the microscopic Hubbard-type 
model to the effective quantum link model, we have con- 
sidered a minimal setup of four lattice sites in the 5=1 
case, as illustrated in Fig. [4^i. We have compared the 
low-energy spectrum of the microscopic and the corre- 
sponding effective Hamiltonian. The spectra coincide for 
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FIG. 4. [color online] a) Pictorial representation of four 
gauge invariant states coupled in second order perturba- 
tion theory in the S = 1 case, b) Expectation values of 
the charge difference ip x i^ x — V^+i^+i an d the electric flux 
E XjX +i as a function of the real-time r, starting from one 
of the four states, under the evolution of the microscopic 
Hamiltonian for t F — U/20 (solid line). The effective quan- 
tum link model (dashed line) yields an expectation value 
— | [l — cos (V2t F tBr/U)] for the product of both opera- 
tors (left plot), while the product of both expectation values 
is — | cos (V2t F t B r/U) [cos (V2t F t B r/U) - l] (right plot), 
signaling the collective dynamics of the coherent fermion hop- 
ping mediated by the quantum link. In the plots t F — 2±b- 

U ^> t Fl and even in the intermediate interaction 
regime U ~ 10 tp the deviations are of order 1%. In view 
of experimental realizations, another relevant question is 
to what extent the Hamiltonian preserves the Gauss law. 
We have considered several initial gauge invariant states, 
evolving in time according to Hjj + AH. For U = 10 tp 
(20 tp), the probability to leave the gauge invariant sub- 
space is below 10% (2%) even for time scales of order 
r ~ 5000 t -1 . Moreover, as demonstrated in Fi^D, the 
expectation values of and E XiX +i display oscillations 
typical of a coherent two-body process, in direct analogy 
with the double well experiments in [TJ [2] . 

At the many-body level, we have also studied the 
flux configuration in the ground state of the microscopic 
Hamiltonian compared to the emergent QLM using ex- 
act diagonalization with open boundary conditions in the 
S = \ representation. As illustrated in Fig. 3a of the 
main text, the microscopic model compares very favor- 
ably with its QLM analog, and gauge invariance is effec- 
tively realized (see Fig. 3b of the main text). 



Implementation procedure 

An experimental implementation will require three ba- 
sic steps: preparation of an initial gauge invariant state, 
evolution via quantum link dynamics, and measurement 
of relevant physical observables. The first step can be 



implemented by preparation of Mott insulator states of 
bosonic and fermionic species on different lattice sites via 
loading in deep optical lattice potentials. Subsequently, 
the ground state or quench dynamics can be realized via 
adiabatic or rapid lowering of the depth of the optical 
lattices. Coherent evolution according to the QLM may 
be validated at the few-body level by performing dou- 
ble well experiments (corresponding to a single quantum 
link) along the lines of Ref. [TJ[2]. A numerical case study 
is presented in the supplementary information. Finally, 
in-situ site-resolved imaging of bosonic particle number 
distributions [3-6 allows one to measure E XyX +i and to 
reconstruct the spin-flux configuration and, thus to quan- 
titatively probe the system. 

Remarks on statistics 

As far as ID implementations are considered, the 
statistics of all basic components does not play a funda- 
mental role (this is not true in 2D and 3D, where Fermi 
statistics is a key-ingredient of matter fields). In particu- 
lar, the gauge bosons may be substituted by fermions (as 
far as S = 1/2 representations are concerned) and mat- 
ter fields may be represented by hard-core bosons. Such 
flexibility may be relevant when considering specific ex- 
perimental implementations, consistently enlarging the 
set of candidate systems. 

Model II: Quantum link models in dipolar systems 

Here we illustrate an alternative route toward realiz- 
ing U(l) quantum link models where a single dipolar 
bosonic species is sufficient to realize a gauge covariant 
link structure. For simplicity, we focus on the S — \ 
setup, although larger spins can, in principle, be achieved 
by considering on-site bosonic interactions. While this 
construction extends to higher dimensions in a straight- 
forward manner, here we discuss the simpler ID imple- 
mentation. 

The microscopic model studied here uses a mixture of 
fermionic and bosonic particles in the presence of strong 
dipolar interactions [TJ [8] . Possible experimental realiza- 
tions are quantum gases of magnetic atoms like Cr [9], 
Er [10], or Dy [nj[l2], and dipolar molecules P3H22]. As 
sketched in Fig. [5j the mixture is confined to a lattice, 
where fermions occupy sites labeled by x and bosons are 
defined on the link sites (x, L) and (x, i?), to the left and 
to the right of x. Bosons can hop only between sites 
(x, R) and (x + 1, L), and serve as natural link variables 
when expressed in terms of Schwinger bosons 

U x ,x+i = b x+1L b x ,R, E x ^ x+ i = - (n x+ i,L — ti x ^r{12) 

Note that here the bosonic index R, L is related to the 
lattice configuration, and, in contrast to Eq. (6) of the 
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main text, it is not associated with an internal degree 
of freedom. By identifying b Xy L = b x: b x ^ R = b x (for x 
even) and b Xy L = b x , b XyR = b x (for x odd), we can relate 
the bosons with spatial indices (x, L) and (x, R) to the 
bosonic species a = 1,2 that arise in model I discussed 
in the main text. The microscopic Hamiltonian of model 
II takes the form 



a) 



L R 



Hdip = -tF ^{^l^x+i + h.c.) + m ^(-1) 



-t B ^2(bl iL b x + liR + h.c.) + u a (-l) x n x , a 

x x,a=L,R 

+ 2U n x,aVa(3[x,y]n y ^ 
x<y a,(3=L,R 

n xi n x,L +n XtR ], (13) 

X 

where v a p[x,y\ is given by the dipolar interaction be- 
tween the particles. Its strength is normalized such that 
vrl[%,% + 1] = 2vlr[%,%] — 1- The Bose-Fermi interac- 
tion is a combined effect of both dipolar and short-range 
potentials. The latter stems from the partial overlap be- 
tween the single-site Wannier functions of bosons and 
fermions in (x,L,R) and x, respectively. In analogy with 
model I, and taking into consideration the fast spatial 
decay of dipolar inter act ions x one can reformulate^ the 
Hamiltonian of Eq. (13) as H dip = U^ X G 2 X + A# dip , 
where 

G x = n% + n XjL + n XjR + X - [(-If + 1] . (14) 

By choosing ujp = U + m, Wfb = ^l,r = U, one then 
obtains a quantum link model with t = 2tB~tF/U, with 
additional gauge invariant diagonal terms generated by 
the dipolar interaction beyond nearest-neighbor sites. 

In case of magnetic atoms [9HTT], the interaction 
regime U > 10 tp may be achieved by properly tun- 
ing the interspecies scattering length and the optical lat- 
tice depth, leading to typical energy scales of the order 
5 nK for, e.g., Dy bosonic gases confined in an optical 
lattice with a lattice spacing of about 200 nm [? ] . Polar 
molecules have large electric dipole moments which can 
be aligned by using external electric fields providing suffi- 
ciently strong constraint energies U [8] when loaded into 
optical lattices [22]. A clear advantage of this setup is 
that it can be straightforwardly adapted to 2D, since, in 
contrast to model I, just one bosonic species is required 
regardless of the dimensionality. Furthermore, dipolar 
interactions are isotropic once the dipole moments are 
aligned perpendicular to the 2D plane [TJ [8] . 

The above realization of gauge fields is based on 
bosonic atoms in a double well potential corresponding to 
the familiar external Josephson effect. Another intrigu- 
ing possibility provided by dipolar interactions makes use 
of an internal Josephson effect [23] , where (instead of two 
wells coupled by a tunnel coupling) one considers atoms 




(U + m) 




FIG. 5. [color online] a) Lattice configuration of model II. 
Fermions hop between sites x and x + 1, while bosons hop 
in double wells on each link, b) Lattice configuration includ- 
ing Hamiltonian parameters, c) Alternative scheme employ- 
ing Raman assisted tunneling rates: fermions with spin are 
loaded into a species-dependent optical lattice, while bosons 
are confined in a periodic potential with half the wavelength. 



with two internal states coupled by a Rabi frequency. 
Such an internal Josephson junction provides an alterna- 
tive realization of quantum link spins. The setup that 
we propose is illustrated in Fig. [5] The gauge invariant 
dynamics are realized by combining a fermionic species 
with spin moving in a spin- dependent optical lattice, and 
a bosonic dipolar species confined in a deep optical lattice 
of half the wavelength. Tunneling of fermions between 
adjacent lattice sites is implemented by a Raman assisted 
transition in the spirit of [24, 25 , coupled to the internal 
Josephson system with the corresponding bosonic link 
site (c.f. Fig. [5]), leading to an effective Hamiltonian term 
of the form ijj[U XiX +itl) x +i. Gauge invariance is then im- 
plemented by considering state-dependent dipole-dipole 
interactions between the bosonic links, and a finite detun- 
ing shift for the Raman transition S. A detailed study of 
this implementation scheme will be reported elsewhere. 
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STRING BREAKING AND FALSE VACUUM 
DECAY IN THE (1 + 1)-D U(l) QUANTUM LINK 
MODEL WITH STAGGERED FERMIONS 

Introduction 

String breaking is an important dynamical mechanism 
in Quantum Chromodynamics (QCD). A quark and an 
anti-quark separated by a distance r are connected by a 
string of color electric flux. The string costs energy in 
proportion to its length, and thus leads to a linearly ris- 
ing quark-anti-quark potential V(r) = ar, where a is the 
string tension. When the quark mass is taken to infin- 
ity, i.e. when one considers an SU (3) Yang-Mills theory 
of pure gluons, in which quarks and anti-quarks appear 
just as external sources of color flux, the confining string 
is unbreakable. In QCD with dynamical quarks, however, 
the string can break by quark-anti-quark pair creation, 
as soon as the energy of the string is larger than the 
mass of the mesons generated in the pair creation pro- 
cess. The dynamics of string breaking, i.e. its evolution 
in real time, is inaccessible to classical simulations us- 
ing lattice gauge theory, which is usually formulated in 
Euclidean time. As for other quantum systems with a 
large number of degrees of freedom, studying the real- 
time evolution of a lattice gauge theory leads to a very 
severe sign problem. For this reason, understanding the 
real-time dynamics of gauge theories from first principles 
remains inaccessible to the traditional methods of theo- 
retical physics. Quantum simulators are a very promising 
tool for overcoming this problem. String breaking is only 
one example of a dynamical mechanism whose real-time 
evolution is interesting but very difficult to investigate. 
The real-time evolution of heavy ion collisions is another 
example, which may serve as a strong motivation and an 
ultimate long-term goal for the development of quantum 
simulators for gauge theories. 

Another interesting phenomenon of great relevance in 
particle physics is the spontaneous breaking and restora- 
tion of symmetries. For example, in the early universe, 
chiral symmetry was intact before it got spontaneously 
broken at sufficiently low temperature. In heavy ion col- 
lisions, one attempts to recreate a droplet of the early 
unbroken phase and observes its evolution back into the 
vacuum state. The time evolution of metastable false 
vacuum states has been discussed in doomsday scenar- 
ios of cosmic evolution, in which the rapid expansion of a 
bubble of true vacuum wipes out the preexisting false vac- 
uum. In one spatial dimension, string breaking by pair 
creation and false vacuum decay are similar phenomena, 
with the dynamical quark and anti-quark playing the role 
of the walls of the true vacuum bubble. In fact, as we will 
see, we can use the same model Hamiltonians to quantum 
simulate both string breaking and false vacuum decay. 

In Wilson's formulation of lattice gauge theory, the 
fundamental gauge degrees of freedom are parallel trans- 



porter matrices taking values in the gauge group, which 
are associated with the links connecting neighboring lat- 
tice sites on a 4-dimensional Euclidean space-time lat- 
tice. In lattice QCD the link variables are 3x3 uni- 
tary matrices with determinant 1, which constitute an 8- 
parameter family of continuously varying classical gauge 
degrees of freedom per link. As a consequence, in Wil- 
son's formulation of lattice gauge theory, there is an 
infinite-dimensional Hilbert space already for each indi- 
vidual link, which largely complicates the construction of 
quantum simulators. 

Quantum link models are an alternative formulation 
of lattice gauge theory, in which the fundamental gauge 
degrees of freedom are represented by discrete quantum 
variables, so-called quantum links, which have a finite- 
dimensional Hilbert space per link. In the simplest case 
of a U(l) gauge theory, a quantum link is simply repre- 
sented by a quantum spin. When one chooses spin |, a 
single qubit per link is hence sufficient to represent the 
gauge field. In an SU (2) gauge theory the quantum link 
is a generalized quantum spin with a continuous gauge 
symmetry, which has at least 4 discrete degrees of free- 
dom. Hence, two qubit s per link are required to represent 
a non-Abelian gauge field. A U (3) quantum link model 
requires 6 and an SU (3) quantum link model requires 20 
discrete states per link, which amounts to 3 or 5 qubits 
per link, respectively. For simplicity, in this supplemen- 
tary material we limit ourselves to U(l) quantum link 
models in one spatial dimension, considering both the 
spin \ representation (one qubit) and the spin 1 repre- 
sentation, with 3 discrete states per link. 

In this supplementary material, we discuss quantum 
link models with a small number of degrees of freedom, 
which still display the physical phenomenon of string 
breaking or false vacuum decay. The proposed model 
Hamiltonians are not yet fully realistic from a particle 
physics point of view. In particular, they address the 
physics of the Schwinger model (i.e. QED in 1 + 1 space- 
time dimensions) rather than QCD in 3 + 1 dimensions. 
For simplicity, the Schwinger model is formulated with 
so-called staggered fermions, which have only one degree 
of freedom per lattice site. It should be pointed out that 
we are not yet addressing the Schwinger model in the con- 
tinuum limit. Hence, several steps will have to be taken 
in order to turn the proposed quantum simulations into 
something that becomes truly useful for particle physics. 
Still, the pathway towards this ultimate goal is clearly 
visible. 



From Classical Background Gauge Fields to 
Quantum Links 

In order to better understand the theoretical frame- 
work of quantum link models, let us begin to discuss the 
familiar situation of fermions moving in a classical back- 
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ground magnetic field described by a vector potential A 
with B = V x A. For simplicity, we consider spinless 
fermions hopping on a lattice with sites x. The fermion 
creation and annihilation operators obey standard anti- 
commutation relations 

{^ y } = 6 xv , {^,^} = {^,^} = 0. (15) 

The term in the Hamiltonian that describes gauge covari- 
ant hopping from a site y to a nearest-neighbor site x is 
then given by iplu xy ip y . Here u xy = exp(i(p xy ) G U(l) is 
a classical gauge parallel transporter associated with the 
link connecting x and y that is obtained by integrating 
the vector potential along the link 



^Pxy 



f 

J x 



dl-A. 



(16) 



After a gauge transformation A' = A — one obtains 

f y ■+ - - 
<Pxy = dl-(A- Va) = (p xy + a x - a y , (17) 

J X 

and hence the parallel transporter transforms as 

u' xy = exp(i(p' xy ) = exp(ia x )u xy exp(-ia y ). (18) 

In order to render the hopping Hamiltonian gauge invari- 
ant, the gauge transformation of the classical background 
field must be accompanied by the gauge transformation 
of the fermion creation and annihilation operators 

i/>' x = exp(ia x )i/; x , ^ x ' = ^ exp(-ia x ), (19) 

which leaves their anti-commutation relations un- 
changed. 

In this work, we are not limiting ourselves to classical 
background fields. Instead the electromagnetic field is 
treated as a dynamical entity with its own quantum dy- 
namics that is intimately connected to the dynamics of 
the fermions. A dynamical lattice gauge field U xy is no 
longer given in terms of a classical background field A. 
Instead, the link variable U xy is an independent quantum 
degree of freedom, whose canonically conjugate "momen- 
tum" variable E xy represents an electric field operator 
that obeys the commutation relation 



\^xy> Uxy] U X y. 



(20) 



In a gauge theory one uses redundant gauge variant vari- 
ables to describe the gauge invariant physics. In order to 
eliminate the unphysical gauge variant states, one must 
impose Gauss' law. Physical states are gauge invari- 
ant, i.e. G x \Sfr) = 0, where 



g x = - £ (/•-• ,, ., - /•:, ,..,) 



(21) 



is the infinitesimal generator of gauge transformations. 
A general (non-infinitesimal) gauge transformation with 
parameters a x is described by the unitary transformation 



V = Y[exp(ia x G x ), 



(22) 



which indeed acts as 



U' 

xy 



V ] U xy V = exp(ia x )U xy exp(-zo^), 
ip f x = V^ip x V = ex.p(ia x )ifj x , 
^=VHlV = ^UM^x)- (23) 



The Hamiltonian of a lattice gauge theory with a dy- 
namical U(l) gauge field coupled to dynamical fermions 
takes the form 



H =~t £ tytUsyt/ty + h.C.) + 9 -Y J El y 



(xy) 



(xy) 



4<? 



\ J2 (U wx U X yU yx U xw + h.c). (24) 



(wxyz) 



Here t is the hopping amplitude and g is the gauge cou- 
pling that determines the electric and the magnetic field 
energies. By construction, the Hamiltonian is gauge in- 
variant, i.e. [H, G x ] — 0. 

In the Hamiltonian formulation of Wilson's U(l) lat- 
tice gauge theory, the link variables are still complex 
phases U xy = exp(ip xy ) G U(l), but the link an- 
gles (p xy are now independent degrees of freedom (un- 
related to a background field A). Correspondingly, the 
canonically conjugate electric field operators are given 
by E xy = —id/d(p xy , such that the commutation rela- 
tion of eq.(20) is indeed satisfied. Since Wilson's parallel 



transporter is a continuous variable, the corresponding 
Hilbert space is infinite-dimensional even for just a sin- 
gle link. In the quantum link formulation of U(l) lattice 
gauge theory, on the other hand, the dimension of the 
link Hilbert space is finite, and given by the 25 + 1 states 
of an integer or half-integer quantum spin S xy on each 
link. In this case, the electric field operator is given by 
E xy = S xy with eigenvalues — 5, . . . , 5, while the gauge 
field is represented by the quantum link operators 

U X y — S xy = S xy -\-iS xy , U^ y = S xy = S xy —iS xy , (25) 

which act as raising and lowering operators of electric 
flux. It is to be noted that the operators E xyi U xyi Ul 



satisfy a local SU(2) algebra at each link (see Eq. 20). 
This allows the realization of a local U(l) gauge invari- 



xy 



ance of the Hamiltonian in Eq. 24 
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Lattice Hamiltonian for the Schwinger Model with 
Staggered Fermions 

The Hamiltonian for a simple quantum link model is 
given by 



H=-t 



(26) 



This model describes a single species of staggered 
fermions (i.e., the fermion field operator ip has only a 
single component) minimally coupled to a U(l) gauge 
field. The electric flux operator defined on each link is 
Ex,x+i = &x x+ii an d the gauge field is represented by 



quantum links U x , x +i = 



°x,x+l 



and U, 



which satisfy the commutation relation 



t 

x,x-\-l 



= s. 



[Ex,X-\-l 1 Uy^y-\-l\ $ XyU x ,X-\-l ' 



x,x-\-l ' 



(27) 



In this supplementary material, we will consider quan- 
tum links both in the spin | and in the spin 1 representa- 
tion. In the minimal spin | representation, the allowed 
electric flux eigenvalues are ±|, i.e. it is impossible to 
have zero flux. This may seem unnatural from the point 
of view of particle physics. Still, as we will see, the result- 
ing physics resembles that of a non-zero vacuum angle 0. 
The spin 1 representation allows the three flux eigenval- 
ues and ±1. 

The factor ( — l) x in the fermion mass term is due to the 
use of staggered fermions. It explicitly breaks the trans- 
lation symmetry by one lattice spacing. This symmetry 
plays the role of a discrete Z(2) chiral symmetry for stag- 
gered fermions. While in the spin 1 case the dynamics of 
the gauge field is governed by the term proportional to 
the gauge coupling g 2 , this term has a trivial contribution 
in the spin \ case. 



Symmetries of the Hamiltonian 

The following symmetry transformations are interest- 
ing to consider: 

• Gauge Invariance: The Hamiltonian commutes 
with the infinitesimal generators of local U(l) 
transformations 



Gx — ^x^x — Ex,x+l + Ex-l,x- 



(28) 



• Chiral Symmetry: The model admits the following 
discrete chiral symmetry transformation: 

X ^x = l/>x+l, = ^l+l , 

X U XyX +l — U x +l,x+2, X E X ,x+l = ^+1,^+2- (29) 

This transformation preserves all terms in the 
Hamiltonian except the mass term. 



• Charge conjugation C: This is again a discrete sym- 
metry, which is implemented as: 

c 4, x = (-ir + vi + i, = (-i) x+ Vx+i, 

C U X)X +i — ^1+1,^+2' C E x , x +i — —E x +i, x +2' (30) 
This transformation leaves all terms invariant. 

• Parity P: The parity transformation is implemented 
as: 

P U x ,x+l — ^- x -\-x^ P E x ,x+l — ~^-x-l,-x- (31) 

• CP Symmetry: The combined symmetry takes the 
form: 



,x+i — U—x,—x+i) Ex,x+i — E—x,—x+i- (32) 



Implementation of Gauss' Law 

To set the conventions, we will use the pictorial repre- 
sentation for the quantum links shown in Figure [6j 



FIG. 6. Left: Spin | quantum link carrying an electric flux 
of +| or — |, respectively. Right: Spin 1 quantum link with 
a flux of +1, 0, or -1, respectively. 

Before we discuss the implementation of the Gauss law, 
note that the staggered fermions introduce a "staggering" 
in the occupation numbers of the vacuum, which resem- 
bles a half-filled state in condensed matter physics, even 
in the absence of the gauge field. Because of the mass 
term, in the vacuum (i.e. the Dirac sea) the fermions oc- 
cupy the odd sites (for m, t > 0, m ^> t). This means that 
the Gauss law cannot be satisfied in the usual way: even 
with equal amounts of electric flux leaving and entering 
a site, the flux divergence G x — ty\$x ~ E x , x +i + Ex-i,x 
depends on whether the site is even or odd. This is an 
unavoidable feature of the staggered fermion formulation. 
We define the Gauss law such that the configuration with 
the staggered occupation in the absence of electric flux 
is a physical state that satisfies it, although not in the 
usual form. All physical states are then required to 
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satisfy 

G X \V) = 



i/jIi/j x - E X)X+ i + E x -i iX + 



("l) a 



=0.(33) 



The vacuum state for spin 1 and t = 0, which obeys 
Gauss' law, is illustrated in Figure [7k. 



String Breaking 

The string formed between an external static charge- 
anti-charge pair is an interesting physical object to be 
studied using quantum simulation. Let us consider string 
breaking in the case of spin 1 quantum links. Figure [7}} 
shows an initial state which has a string joining an exter- 
nal static charge Q and an anti-charge Q placed at both 
ends. As the fermions become lighter, the sequence of fig- 
ures ^jp - [7^ shows the breaking of the string due to the 
creation of fermion-anti-fermion pairs, and the evolution 
to the final state of two mesons at the ends. For the static 
t = case, the vacuum energy is given by Eo = — mL/2, 
while the energies of the initial unbroken string state and 
the final two-meson state take the form 



string 



-(L-l), 



^mesons ~ E = 2( — + Tfl) . 



(34) 



The critical distance for string breaking is determined by 



the condition E< 



string " 



-E, 



results in L = Am/g 2 



mesons 

3. 



- {L — 3) — 2m and thus 



False Vacuum Decay in a Minimal Model with C 
and P Symmetry 

The minimal model with S = \ quantum links may 
seem unnatural from a particle physics point of view, 
because it does not allow zero flux. However, it resembles 
a non-zero vacuum angle 6 = tt. In order to explore these 
dynamics, let us consider the Hamiltonian 

X 

+m^(-lf4^-2a^(-lf^ +1 . (35) 
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FIG. 7. a) vacuum, b) string induced by a static external QQ 
pair, c) broken string, d) evolution, e) final state with two 
mesons separated by vacuum. 



on the other hand, charge conjugation and parity remain 
exact symmetries, as it is the case here. Interestingly, 
in the Schwinger model at 9 = 7r, both charge conjuga- 
tion and parity are spontaneously broken. As we will 
see, this is also the case in the model that we have just 
introduced, at least for sufficiently large values of the 
mass m (and for t^O), while for small m, charge conju- 
gation and parity are restored. Hence, the model with 
the staggered flux term proportional to a mimics the 
physics of = 7r, which is again very difficult to sim- 
ulate with a classical computer, due to a very severe sign 
problem. Let us finally notice that, after implementing 
Gauss law, the staggered electric field term can be recast 
into a renormalization for the fermionic mass term, that 
is m e ff = m — cr, so that one recovers exactly the model 
in Eq. (3) of the main text. 

Let us consider the physics in the vacuum sector, using 
periodic boundary conditions of a finite system with an 
even number of lattice sites L. For t = 0, there are three 
candidate vacuum states, which are depicted in Figures 
[8] as well as and (9}d. First of all, it is instructive to 



The last term favors a staggered flux pattern and is 
charge conjugation and parity invariant. It does, how- 
ever, explicitly break chiral symmetry. Since chiral sym- 
metry is anyways explicitly broken by the Gauss law im- 
plementation, this additional breaking is unimportant. 
Terms linear in the electric field are characteristic of a 
non-zero vacuum angle 0. For ^ 0,7r, charge conjuga- 
tion and parity are both explicitly violated. At = 0, 7r, 



FIG. 8. C and P invariant candidate vacuum state of the spin 
| model. 

convince oneself that these three states indeed satisfy the 
Gauss law. The state shown in Figure|8]is invariant under 
both charge conjugation and parity, while the two states 
in Figure [9] are charge conjugation and parity partners of 
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= 7r. It would be interesting to quantum simulate the 
decay of the false vacuum of Figure [8] into the true vacua 
of Figure [9] in real time. The quench dynamics discussed 
in the main text is directly related to this phenomenon. 



FIG. 9. Competing vacua of the spin 
and P partners of each other. 



model, which are C 



each other. At t 
given by 

£ = -°L ■ 



the energies of the three states are 



m-, E' = E% = 



-m 



2' 



(36) 



such that E f — E = Eq — E = (a — m)L. Hence, for 
m < a (with a > 0) the state depicted in Figure [8] is the 
true vacuum (i.e. the Hamiltonian eigenstate of lowest 
energy respecting the Gauss law). In this case, both C 
and P are unbroken. On the other hand, for m > cr, the 
two states shown in Figure [9] have a lower energy and 
thus C and P are then spontaneously broken. This is 
reminiscent of the situation in the Schwinger model at 



Conclusions 

As we have seen, a simple (l+l)-d quantum link 
model with staggered fermions allows us to address in- 
teresting dynamical questions using quantum simulation. 
Along the way towards addressing similar problems in 
QCD, one faces the challenges of higher dimensions, non- 
Abelian gauge fields, and multi-component fermions. In 
the quantum link formulation of QCD all these elements 
are present, and we see no fundamental obstacle against 
implementing then with ultra cold matter in optical lat- 
tices. As one goes beyond (1+1) dimensions, gauge field 
acquire transverse degrees of freedom and display a richer 
dynamics. Since quantum link models are gauge invari- 
ant by construction, these dynamics are those of QCD in 
the continuum limit. How to best realize this in quan- 
tum simulations will be addressed in forthcoming publi- 
cations. 



